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ABSTRACT 

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed 
with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO (2-1) emission of four 
massive cores. We detect dust condensations down to ~ 0.3 M and separate millimeter sources down to 0.4" or 
< 1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The 
CO (2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known 
outflow directions. This, in combination with additional cores from the literature observed at similar mass 
sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 
3 orders of magnitude. Among the 18 regions, ~ 30% show no signs of fragmentation, while 50% split 
up into > 4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as 
bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with 
the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and 
turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic 
simulations, and obtained that the low-fragmented regions are well reproduced in the magnetized core case, 
while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic 
field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the 
initial magnetic field/turbulence balance in each particular core. 

Subject headings: stars: formation — ISM: individual objects (IRAS 22172+5549, IRAS 22134+5834, 
IRAS 22198+6336, AFGL5142) — radio continuum: ISM 



1. INTRODUCTION 

A well-observed trend is that most stars form in clus- 
ters/groups, rather than in isolation (e. g., Carpenter 2000; 
Lada & Lada 2003), and thus the understanding of cluster 
formation is intimately related to the understanding of star 
formation. One of the results found from statistical studies 
of clusters at optical/infrared wavelengths, carried out down 
to mass sensitivities of 0.1-0.5 M and spatial resolutions of 
~ 1000 AU, is that there is a correlation between the mass of 
the most massive star in the cluster and the number of mem- 
bers of the cluster (e. g., Hillenbrand 1995; Testi et al. 1999; 
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Qiu et al. 2008). However, it is not clear which is the physical 
reason for this correlation: it could either be the result of the 
initial mass function (e. g., Testi et al. 2000; Parker & Good- 
win 2007), or it could indicate that the most massive stars 
need a rich cluster to be formed (e. g., Stahler et al. 2000; 
Krumholz & Bonnell 2007). Alternatively, the most massive 
stars may play an important role to increase the number of 
members in their surroundings, by inducing the formation of 
new stars (e. g., Hester & Desch 2005; Deharveng et al. 2012). 
Thus, it is necessary to investigate whether the mentioned cor- 
relation is already established at the beginning of the forma- 
tion of the cluster, or it is a result of cluster evolution. 

To answer these questions about cluster formation it is of 
fundamental importance to characterize clusters in their very 
first stages of formation, when they are deeply embedded 
in a massive dense core and thus detectable in the millime- 
ter/submillimeter range. However, in general the present mil- 
limeter/submillimeter instrumentation does not allow to reach 
the required mass sensitivity and spatial resolution (< 0.5 M , 
~ 1000 AU), given the typical distances of massive dense 
cores of > 1 kpc, in a reasonable amount of time for a sample 
of regions. Thus, several studies have focused on millime- 
ter observations of individual massive cores, associated with 
intermediate/high-mass protostars, down to spatial scales of 
< 1000 AU and reaching mass sensitivities around 0.1 M Q 
(e.g., Beuther et al. 2004; Zapata et al. 2005; Fuente et 
al. 2005; Leurini et al. 2007; Neri et al. 2007; Teixeira et 
al. 2007; Rodon et al. 2008; Su et al. 2009; Tang et al. 
2010), but no studies up to date focus on the formation of 
clusters from a statistical point of view. More recently, mas- 
sive cores at even earlier evolutionary stages (being dark or 
'quiet' in the infrared -the so-called massive infrared dark 
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clouds), have been imaged at high angular resolution to in- 
vestigate fragmentation (e. g., Swift 2009; Zhang et al. 2009; 
Bontemps et al. 2010; Pillai et al. 201 1; Zhang & Wang 201 1; 
Wang et al. 2011; Wang et al. 2012b). In particular, Bon- 
temps et al. (2010) imaged with the Plateau de Bure Inter- 
ferometer (PdBI) in compact (CD) configuration six IR-quiet 
massive cores in Cygnus X, studying fragmentation down to 
1400 AU, and find that the fragmentation level is higher than 
expected for a turbulence regulated monolithic collapse sce- 
nario (e.g., McKee & Tan 2002, 2003), but lower than ex- 
pected for a pure gravo-turbulent scenario (e. g., Klessen et 
al. 2005; Bonnell et al. 2007; see Csengeri et al. 2011 for 
further details on both scenarios). This calls for further stud- 
ies to understand the fragmentation process. In this paper we 
report on observations carried out with the PdBI operated in 
A configuration toward four massive dense cores associated 
with intermediate/high-mass protostars. These PdBI millime- 
ter observations provide a spatial resolution and sensitivity 
comparable to those reached in the optical/infrared studies. 
These four regions, completed with the regions from the liter- 
ature mentioned above, constitute a step towards a statistical 
study of the first stages of cluster formation. 

The four observed regions were selected from the litera- 
ture following these criteria: i) regions must have been ob- 
served with a millimeter/sub millimeter interferometer in a 
compact configuration (with angular resolution < 7"), and 
must show a compact condensation of gas and/or dust, to as- 
sure that the intermediate/high-mass protostar is deeply em- 
bedded in a massive dense core; ii) regions must be at a dis- 
tance < 3 kpc, to obtain the required high spatial resolution 
and high mass sensitivity; iii) regions must span a range of 
luminosities between ~ 300 L Q and ~ 10 s L , which is the 
range in which clustering becomes dominant in the star for- 
mation process; iv) regions must belong to different cloud 
complexes to avoid biases due to similar initial conditions; 
and v) regions had to be observable from the PdBI when the 
A configuration was operative. These selection criteria will 
allow us to study the richness of the protoclusters as func- 
tion of the luminosity and compare it to the optical/infrared 
results found for more evolved clusters. The four selected re- 
gions are IRAS 22198+6336 (122198); a region ~ 15" north 
of IRAS 22172+5549 (I22172N); AFGL5142 (A5142); and 
IRAS 22134+5834 (122134). 122198 is a Class source of 
~ 300 L Q harboring a hot core and driving a quadrupolar out- 
flow (Sanchez-Monge et al. 2010). I22172N and A5142 have 
luminosities of ~ 800 and 2200 L©, respectively, are also clas- 
sified as Class sources, and are associated with outflows 
(e.g., Hunter et al. 1999; Fontani et al. 2004; Zhang et al. 
2007). Finally, 122134, with a luminosity of ~ 12000 L , 
is an UCHII region embedded within a dust condensation 
(Sanchez-Monge 201 1). Thus, all the selected regions harbor 
intermediate/high-mass protostars deeply embedded in mas- 
sive dense cores, being excellent protocluster candidates. We 
will follow the nomenclature of Williams et al. (2000) and 
Bontemps et al. (2010) where a massive dense core refers to a 
dense gas structure of ~ 0.1 pc of size and > 20 M of mass 
which does not necessarily collapse into one star, but can frag- 
ment in compact condensations and form a small cluster of 
stars. 

In Section 2 we describe the PdBI observations in the most 
extended (A) configuratiorQ in Section 3, we present the re- 

12 The PdBI A-configuration dataset of 1.3 mm continuum and CO (2-1) 
emission presented in this work allowed us to simultaneously detect emission 



TABLE 1 

Parameters of PdBI 1.3 mm continuum observations 





Phase center 


synth. beam 


P.A. 




Source 


a(J2000) 


<5(J2000) 


("X") 
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rms a 


I22172N 


22:19:08.60 


+56:05:02.0 


0.47 x 0.38 


26 


0.48 


122134 


22:15:09.23 


+58:49:08.9 


0.56 x 0.49 


41 


0.33 


122198 


22:21:26.78 


+63:51:37.6 


0.43 x 0.39 


27 


2.0 


A5142 


05:30:48.02 


+33:47:54.5 


0.52 x 0.36 


12 


2.8 



Rms noise in mjy beam 



suits of the 1.3 mm continuum and CO (2-1) emission for the 
four regions; in Section 4, we present the total sample of 18 
massive dense cores (14 from literature + 4 from this work); in 
Section 5 we discuss the role of the global properties of mas- 
sive dense cores in the fragmentation process, and compare 
our observations to radiation-magneto-hydrodynamic simula- 
tions of collapsing massive dense cores. Finally, in Section 6 
we summarize our main conclusions. 

2. OBSERVATIONS 

The PdBI observations, partially described in Palau et al. 
(201 1), were carried out on 2010 January 10th with the array 
in A configuration at 230.538 GHz to observe A5142. The 
baselines range from 136 to 760 m. The system temperature 
was 150-200 K and the averaged atmospheric precipitable 
water vapor was ~ 2 mm. Radio frequency bandpass cali- 
bration was carried out using B085 1+202. Gain phases and 
amplitudes were calibrated by observing the nearby quasars 
B0552+398 and J05 12+294, yielding a phase rms of 20-60°. 
The adopted flux density for the flux calibrator, 3C273, was 
16.57 Jy. A subsequent track was carried out on 2010 January 
18th to observe in track-sharing mode the I22172N, 122134, 
and 122198 regions. The system temperature was around 
150 K and the precipitable water vapor was ~ 1 mm. 3C454.3 
was used for the radio frequency bandpass calibration. Phases 
were calibrated with B2146+608, and B2037+511, and we 
obtained a phase rms of 20-40°. A flux of 0.23 Jy was fixed 
for the flux calibrator, B2146+608. The absolute flux den- 
sity scale was determined from 3C273 and B2 146+608, with 
an estimated uncertainty around 30%. The absolute position 
accuracy was estimated to be around ~ 0.05". A correla- 
tor unit of 40 MHz of bandwidth with 256 spectral channels 
was used to observe the CO (2-1) line in each polarization. 
Three additional units (per polarization) of 320 MHz band- 
width with 128 channels were used to observe the continuum 
across ~ 1 GHz. 

Calibration and imaging were performed using the Cue and 
MappinJI packages of Gildas software, using the standard 
procedures. The continuum wu-data of 122134 and the CO(2- 
1) Mu-data of all the four sources were tapered to 400 m base- 
lines (Mi>-distance where the weight of the taper Gaussian is 
37%), to increase the signal-to-noise (S/N) ratio. For the 
CO (2-1), we smoothed the channel maps to a velocity res- 
olution of 0.6 km s per channel. Typical rms noises (in line- 
free channels) per channel are 20-30 mjy beam 1 . The final 

of complex organic molecules in 122198 and A5142. This is reported in Palau 
etal. (2011). 

13 http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. — PdBI maps at 1,3 mm in a field of view of 20" (~primary beam) of diameter. Top-left: I22172N. Contours are -4, 4, and 8 times the rms noise of 
the map, 0.55 mJy beam -1 . Synthesized beam, shown in the bottom right corner, is 0.47" X 0.38", PA~ 26.0°. Top-right: 122134. Contours are -4, 4, and 8 
times the rms noise of the map, 0.33 mJy beam -1 . Synthesized beam is 0.56" X 0.49", PA~ 40.8°. Green contours correspond to the 3.6 cm emission tracing 
an UCHII region, and are -4, 4, 8, and 16 times 30 /iJy beam -1 (Sanchez-Monge 2011). Bottom-left: 122198. Contours are -4, 4, 8, 16, and 32 times the rms 
noise of the map, 2.0 mJy beam -1 . Synthesized beam is 0.43" X 0.39", PA~ 27.1°. Bottom-right: A5142. Contours are -4, 4, 8, 12, 16, 20, and 24 times the 
rms noise of the map, 2.8 mJy beam -1 . Noise is limited by dynamic range and phase noise. The synthesized beam is 0.52" X 0.36", PA~ 12.0°. In all panels the 
red contour corresponds to the identification level of 6cr. 



continuum maps were built by carefully selecting the chan- 
nels of the 320 MHz units which were free of line emission. 
The phase centers, synthesized beam sizes, and rms of the fi- 
nal cleaned continuum maps are listed in Table Q] 122198 and 
A5142 images are limited by dynamic range and A5142 ad- 
ditionally has phase noise. The primary beam of the PdBI at 
1.3 mm is 22", allowing us to cover a field of view of 0.1- 
0.2 pc. 

3. RESULTS 

3.1. 1 .3 mm continuum 

In Fig. Q] we present the cleaned images of the four regions 
observed with the PdBI in A configuration (hereafter, PdBI- 
A). Sources above 6<x (with cr the rms noise level) were iden- 
tified in the final cleaned images, and are listed in Table [2] In 
order to confirm that the 6cr identifications in the cleaned im- 
ages are well detected also in the wu-plane, we fitted in the uv- 



plane a Gaussian or a point-like source to the strongest source 
identified in the cleaned map, and subtract the fitted source 
to the observed visibilities (wy-plane). Then, we fitted in the 
residual visibilities the second strongest source identified in 
the cleaned map, etc. All the identified sources (in the cleaned 
map) were detected in the t/y-plane with a S/N ratio larger 
than 10 (except for 122134, for which S/N in the MD-plane 
was > 4.6), giving us confidence on the identification method. 
Position, flux densities and sizes for each source were deter- 
mined from the fits of a Gaussian or a point-like source in the 
Mu-plane. These parameters are listed in Table[2] together with 
an estimate of the mass for each millimeter source. 

In the I22172N field we detected four sources: three faint 
point sources and a strong extended source, MM2, which is 
elongated roughly in the north-south direction, with a size of 
2300 x 1000 AU, and a mass of 2-10 M Q (Table 0. The 
mass is derived assuming a dust temperature of 10-30 K, and 
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TABLE 2 

Parameters of the sources detected with the PdBI at 1 .3 mm (continuum) for the 4 regions studied in this work 





Position" 


Deconv.ang.size 


P A 


Deconv. size 


rpeak a 
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c ;i 
J v 


IViaSS 






/ to f\r\rx\ 

o-(J2000) 


c/ to f\r\rx\ 

<5(J2000) 




V > 


TAT TV ATT^ 


\ iiij y ucd.111 j 


(m Tvl 


(M \ 


A cc r\r* i n 1 1 nn 
r\ a S (JCl a LHJ11 


I22172N-MM1 


22:19:08.155 


56:05:11.41 


- 


- 


< 1000 


9.0 ±0.9 


9±4 


0.7-3 


IR 


I22172N-MM2 


22:19:08.974 


56:05:02.97 


0.94 x 0.40 


2 


2260 x 960 


6.1 +0.5 


29 ±9 


2-10 


Outf 


I22172N-MM3 


22:19:09.178 


56:05:02.37 


- 


- 


< 1000 


8.1+0.6 


8±2 


0.6-3 


- 


I22172N-MM4 


22:19:09.199 


56:05:04.66 


- 


- 


< 1000 


6.3 + 0.6 


6±2 


0.5-2 


- 


I22134-MM1 


22:15:08.778 


58:49:07.90 


- 


- 


< 1300 


2.8 ±0.4 


2.8 ± 1.0 


0.2-1 


- 


I22134-MM2 


22:15:09.395 


58:49:13.04 


— 


— 


< 1300 


2.1 ±0.4 


2.1 ±0.8 


0.2-0.8 


_ 


I22134-MM3 


22:15:09.455 


58:49:06.45 






< 1300 


2.2 ±0.4 


2.2 ± 0.7 


0.2-0.9 




I22134-MM4 


22:15:09.804 


58:49:11.33 






< 1300 


2.1 ±0.5 


2.1 ±0.8 


0.2-0.8 




I22198-MM2 


22:21:26.753 


63:51:37.90 


0.68 x 0.39 


-35 


520 x 300 


91.9 ±2.0 


250 ± 70 


0.5-2 


HC,Outf 


I22198-MM2-S 


22:21:26.807 


63:51:37.14 






< 300 


18.8 ±0.8 


18 ±5 


0.1-0.6 


IR,Outf? 


A5142-MM1 


05:30:48.031 


33:47:54.64 


0.68x0.51 


-86 


1200 x 900 


38.3 ±2.8 


210 ±60 


2-9 


HC,Outf 


A5142-MM2 


05:30:48.024 


33:47:53.67 


0.55 x 0.25 


18 


1000 x 400 


62.6 ± 2.8 


150 ±40 


2-6 


HC,Outf 


A5142-MM3 


05:30:47.953 


33:47:52.57 






< 700 


22.5 ± 1.3 


23 ±7 


1^1 


OA,W 


A5142-MM6 


05:30:47.749 


33:47:53.03 






< 700 


21.9 ± 1.3 


22 ±6 


0.9-4 


OA,W 


A5142-MM7 


05:30:47.980 


33:47:55.53 






< 700 


18.6 ± 1.3 


19 ±6 


0.8^t 


IR? 


A5142-MM8 


05:30:48.152 


33:47:57.78 






< 700 


18.6 ± 1.3 


19 ±5 


0.8-4 




A5142-MM9 


05:30:48.257 


33:47:51.52 






< 700 


19.4 ± 1.3 


20 ±6 


0.9^1 


OA 



11 Position and flux density derived by fitting a point source or a Gaussian in the iw-plane. Peak intensity measured directly in the cleaned 
image. Uncertainty in the peak intensity, cr, is the rms noise of the cleaned image for the sources fitted with a Gaussian, or the rms given 
by the fit of a point source. Error in flux density has been calculated as -\j(cr mmcc /6 bmlm ) 2 + (cr flux _ sc . llc ) 2 (Beltran et al. 2001), where 
^source and #beam are the size of the source and the beam respectively, and cr flux _ sca i c is the error in the flux scale, which takes into account 
the uncertainty on the calibration applied to the flux density of the source (S v x % unC enamty)- Peak intensities and flux densities (with 
uncertainties) are corrected for the primary beam response. 

b Masses derived assuming a dust temperature of 10-30 K (except for I22198-MM2, A5142-MM1 and A5142-MM2, for which we used 
a range of 30-100 K, Palau et al. 201 1), and a dust (+gas) mass opacity coefficient at 1.3 mm of 0.00899 cm 2 g~' (Ossenkopf & Henning 
1994, see main text). The uncertainty in the masses due to the opacity law is estimated to be a factor of 2. 

c Association with signposts of stellar activity: IR = infrared source; Outf = source driving outflow; OA = source along outflow axis; HC 
= hot molecular core (Palau et a l. 201 1); W = water maser emission at velocities about 5 km s~' away from systemic velocities (Hunter 
et al. 1995; 1999); see Sect.[3]2]for further details. 



a dust mass opacity coefficient at 1.3 mm of 0.899 cm 2 g _1 
(column 6 of Table 1 of Ossenkopf & Henning 1994, corre- 
sponding to agglomerated dust grains with thin ice mantles 
at densities ~ 10 6 cm -3 ). We estimated an uncertainty in the 
masses due to uncertainty in the dust opacity of about a fac- 
tor of 2. MM3 and MM4 fall within 3" to the (north)east 
(< 7200 AU) of MM2, and are coinciding with a NH3 emis- 
sion peak (for which a rotational temperature of ~ 20 K is es- 
timated, Sanchez-Monge et al., in prep.), while MM1 is ~ 12" 
to the northwest of MM2. For MM1, MM3 and MM4, we es- 
timated masses around 0.5-3 M , making the same assump- 
tions given above. 

Towards 122134 we detected four point sources with a flux 
of ~ 2 mJy, located in all directions and at similar distances, 
4" or ~ 10000 AU, from the UCHII region. The masses of the 
fragments are ~ 0.7 M , and the UCHII region is not detected 
at 1 mm. Such a non-detection is consistent with the extrap- 
olation of the flux density of the centimeter source at 1.3 cm 
(~ 0.8 mJy beam -1 ) to 1.3 mm, using the measured spectral 
index of -0.18 (Sanchez-Monge 201 1). 

For 122198 we detected one strong source, MM2 (following 
the labeling in Sanchez-Monge et al. 2010), extended in the 
southeast-northwest direction, with a faint extension at 6<x, 
MM2-S, 0.8" (600 AU) to the southeast, already reported in 



Palau et al. (2011). MM2 has a size of 500 x 300 AU at a 
PA=-35°, and its flux density, of ~ 250 mJy, corresponds to 
a mass of ~ 0.5-2 M . Here we assumed a dust temperature 
range of 30-100 K because MM2 is associated with emission 
from complex organic molecules tracing a hot core at ~ 100 K 
(Palau et al. 2011). The faint extension MM2-S has a flux 
density corresponding to a mass of ~ 0.5 M (assuming a 
dust temperature of 10-30 K). 

Finally, for the case of A5142, the millimeter emission 
is dominated by two partially extended and strong sources, 
MM1 and MM2, separated by ~ 1" or 1800 AU, which are 
surrounded by five faint point-like sources (labeled following 
Zhang et al. 2007). The masses of MM1 and MM2 are 2- 
9 M (assuming, as in the case of I22198-MM2, a dust tem- 
perature of 30-100 K, Palau et al. 201 1), and both sources are 
associated with water masers. MM1 is elongated in the east- 
west direction, with a size of 1200 x 900 AU, while MM2 is 
oriented in the north-south direction and is sligthly smaller, of 
1000 x 400 AU. Concerning the other sources, they are com- 
pact and faint, with masses estimated around 1-4 M ; MM7 
lies at about 1 " to the northwest of MM 1 , and MM3 and MM6 
are associated with water maser emission (Hunter et al. 1995, 
1999). 
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3.2. CO ( 2-1 ) emission, 2MASS and Spitzer/IRAC sources 

In Figs. [2] [3] @] and [5] we present the 1.3 mm continuum 
emission as in Fig. Q] with the high-velocity blueshifted and 
redshifted CO (2-1) emission superimposed. The CO (2-1) 
emission is heavily filtered out (the PdBI operating at 1 mm in 
A configuration is not sensitive to structures larger than > 1 "), 
and reveals only the high-velocity most compact knots. As 
an example, we estimated that the PdBI-A misses about 60- 
70% of the high velocity CO (2-1) emission in 122198 (see 
below) compared to the flux measured with the SMA in com- 
pact configuration (Sanchez-Monge et al. 2010). Although 
suffering from strong filtering, most of the high-velocity knots 
recovered with the PdBI-A are consistent with the previously 
known outflows in each region, and for I22172N the CO (2-1) 
emission reveals one new outflow. 

In order to further understand the nature of the 1 mm con- 
tinuum sources, we searched for 2MASS and Spitzer/IRAC 
sources (at 3.6 jum) within the 22" PdBI field of view, and the 
positions are marked in Figs. |2] [3] IH and [5] The positions 
of the Spitzer/IRAC sources were extracted from the image 
at 3.6 fim using the Source-Extractor software (Bertin & 
Arnouts 1996). In the last column of Table |2] we list which 
millimeter sources are associated with infrared emission, wa- 
ter maser emission (Hunter et al. 1995, 1999), or hot molec- 
ular cores (Palau et al. 2011). We also indicate if the sources 
are driving outflows or are located along outflow axes. Those 
millimeter sources with no evidences of stellar activity could 
either be starless or protostellar. We favor the protostellar 
option because starless cores have flat column density pro- 
files with typical (half -power) radius of about 3000-6000 AU 
(e.g., Motte et al. 1998; Kandori et al. 2005), implying that 
the PdBI-A should recover a fraction of correlated flux for 
these cores of only < 0.1% (following Palau et al. 2010). 
Thus, starless cores are presumably filtered out in our obser- 
vations. In the following, we give a small description of the 
CO (2-1) PdBI-A emission for each source, together with a 
brief interpretation of the nature of each millimeter source. 

I22172N: — The first interferometric observations studying 
the 3 mm and CO (1-0) emission in I22172N were carried out 
with OVRO by Fontani et al. (2004), who find one millimeter 
condensation with one outflow displaced about 4" to the east. 
Our PdBI-A observations reveal one strong and elongated 
source, MM2, coincident with the peak of emission seen at 
3 mm by Fontani et al. (2004), and which seems to be driving 
a small outflow not reported before (Fig. |2}. The blueshifted 
lobe of this small outflow is associated with a faint H2 feature 
reported in figure 15-bottom of Fontani et al. (2004), which is 
elongated also in the southeast-northwest direction, confirm- 
ing the presence of this outflow. MM2 and its nearby sources 
MM3 and MM4 are deeply embedded in a dense core traced 
in NH3(1,1) (Sanchez-Monge et al., in prep.) and lie in a re- 
gion with no near-infrared emission (see figures 10 and 15 of 
Fontani et al. 2004). On the other hand, the large-scale CO (1- 
0) outflow found by Fontani et al. (2004) seems to be driven 
by a near-infrared and Spitzer/IRAC source, labeled IRS 1 in 
the figure, which falls exactly at the center of symmetry of 
the two CO (1-0) outflow lobes but which has no dust asso- 
ciated. Concerning MM1, it falls 1" to the northwest of one 
bright near-infrared and Spitzer/IRAC source. Since the po- 
sitional uncertainty of Spitzer/IRAC sources is also ~ 1", the 
infrared source and MM1 are possibly associated. The faint 
detection of MM1 with the PdBI-A, together with its strong 
intensity in the near-infrared, and the lack of a strong out- 
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Fig. 2. — I22172N: PdBI-A map of the continuum at 1.3 mm (black con- 
tours as in Fig.fTJ with the CO (2-1) high- velocity blueshifted (blue contours, 
integrated from -52.4 to -44.6 km s - ') and redshifted (red contours, inte- 
grated from -40.4 to -36.8 km s _1 ) emission overplotted (systemic velocity 
is -43.8 km s ). Red and blue arrows mark the directions of the outflows 
seen with the PdBI-A. Grey (dashed, solid) contours and black arrows indi- 
cate the CO (1-0) outflow previously reported by Fontani et al. (2004, with 
a beam ~ 6"). Contours of the PdBI CO (2-1) emission start at 55% and 
increase in steps of 15% of the maximum, 2.13 Jy beam -1 km s~', for the 
blueshifted emission, and start at 60% increasing in steps of 20% of the max- 
imum, 0.87 Jy beam -1 km s~', for the redshifted emission. Plus signs cor- 
respond to 2MASS sources, and crosses to Spitzer/IRAC sources at 3.6 fim. 
The synthesized beam of the PdBI CO (2-1 ) emission, of 0.57" X 0.49" with 
P.A.= 39.2°, is shown in the bottom-right corner. 

flow suggests that MM1 is a slightly more evolved object than 
IRS1 and MM2 (MM3 and MM4), as expected if star for- 
mation in I22172N is triggered by the WR star HD 21 1853, 
which is about 5' to the (north)west (Liu et al. 2012). Overall, 
I22172N harbors a rich variety of young stellar objects: two 
protostellar candidates detected only in the millimeter range 
(MM3 and MM4), one confirmed protostellar object with no 
infrared emission and driving a small outflow (MM2), one 
protostellar object detected only in the infrared and driving a 
larger outflow (IRS1), and one protostellar object detected in 
the millimeter and infrared, with no outflow (MM1). 

122134: — The NH3(1,1) and (2,2) emission of this region is 
studied in Busquet (2010) and Sanchez-Monge (2011), who 
find that the dense gas is mainly concentrated towards the 
southeast of the UCHII region, with faint NH3 emission form- 
ing a tilted 'U' structure around the UCHII region (Fig. 0J, 
where the broadest linewidths and highest rotational temper- 
atures are found. Thus, the dense gas seems to be tracing 
the walls of a cavity which could have been created by the 
UCHII region, similar to the cavity recently found by Fontani 
et al. (2012) in IRAS 20343+4129. Interestingly, the mil- 
limeter sources found with the PdBI-A are distributed along 
the 'U' NH3 structure, at the edges facing the UCHII region, 
suggesting that they could be dust condensations affected by 
the expanding ionization front creating the cavity. As for 
the CO (2-1) emission, all the previous observations are done 
with single-dish and mainly trace the large-scale CO structure 
(e.g., Beuther et al. 2002b), which cannot be related to the 
knots found with the PdBI-A. We find one main redshifted 
knot, in between MM1 and MM3, and one strong blueshifted 
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Fig. 3. — Idem as Fig. [2] for 122134. Grey contours correspond to the in- 
tegrated emission of NFFjd.l) (Busquet 2010, Sanchez-Monge 2011). Con- 
tours of the PdBI CO (2-1) emission start at 45% and increase in steps of 13% 
of the maximum, 1.98 and 2.16 Jy beam - ' km s , for the blueshifted and 
redshifted emission, respectively. CO (2-1) blueshifted emission has been 
integrated from -31.7 to —22.1 km s , while CO (2-1) redshifted emis- 
sion has been integrated from -14.9 to -4.7 km s -1 (systemic velocity is 
-18.3 km s ). The star indicates the position of the UCHII region (Fig. [T] 
and Sanchez-Monge 2011). The synthesized beam of the PdBI CO (2-1) 
emission, of 0.55" X 0.49" with P.A.= 42.5°, is shown in the bottom-right 
corner. 
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Fig. 4. — Same as Fig. |2]for 122198. Grey (dashed, solid) contours cor- 
respond to (red-, blue-shifted) HCN(l-O) outflow emission; and the grey 
dotted contours correspond to the redshifted HCO + (1-0) outflow emis- 
sion (beam around 3", Sanchez-Monge 2011). Contours of the PdBI 
CO (2-1) emission start at 18% and increase in steps of 18% of the max- 
imum, 1.75(1.77) Jy beam - km s , for the blueshifted(redshifted) emis- 
sion. CO (2—1) blueshifted emission has been integrated from -31.0 to 
-19.7 km s" 1 , while CO (2-1) redshifted emission has been integrated from 
-5.3 to 7.9 km s" 1 (systemic velocity is —10.5 km s ). The star marks the 
position of a hot molecular core (Sanchez-Monge et al. 2010; Palau et al. 
201 1), and the cross indicates the position of the 3.6 /im Spitzer/IRAC source 
found in the field. The synthesized beam of the PdBI CO (2-1) emission, of 
0.55" X 0.50" with P.A.= 42.5°, is shown in the bottom-right comer. 



knot about 3" the south. Only 075 to the south of the red- 
shifted CO knot there is a 2MASS and Spitzer/IRAC source, 
as well as emission at Act detected with PdBI-A at 1 mm. All 
this suggests that this infrared source, labeled IRS1 in Fig.[3j 
could be the driving source of an outflow elongated in the 
north-south direction. Observations in the CD configuration, 
recovering larger angular structures, are required to confirm 
the existence of such an outflow. Thus, in this region we found 
4 millimeter sources with no infrared emission plus one faint 
millimeter source being bright in the infrared, with the UCHII 
region located at the center of all these 5 young stellar objects. 

122198:— The PdBI-A CO (2-1) emission in 122198 is re- 
solved into several knots which are mainly aligned along 
the (south)east-(north)west direction, with additional faint 
redshifted knots to the northeast of MM2 and one strong 
blueshifted knot about 5" to the southwest of MM2 (Fig. El). 
This could be interpreted as one single outflow with a wide 
opening angle, or as two independent outflows, outflow A' 
and 'B' as marked in Fig. [4] Observations of the high-velocity 
emission of different outflow tracers (carried out with the 
PdBI in CD configuration, Sanchez-Monge 201 1) seem to fa- 
vor the two-outflow scenario. As shown in grey contours in 
Fig.|U the HCN (1-0) high-velocity emission is bipolar, cen- 
tered on MM2, and aligned in the southwest-northeast direc- 
tion. Similarly to HCN, the SiO (2-1) high-velocity emission 
is only found elongated to the southwest of MM2 (see fig- 
ure 1 of Sanchez-Monge et al. 2010). On the other hand, the 
redshifted emission of HCO + (1-0) is elongated in a direc- 
tion perpendicular to the HCN outflow (see dotted grey line 
in Fig.Q. Since SiO and HCN outflows are more common in 
Class protostars than in more evolved Class I objects (e. g., 
Tafalla et al. 2000; Shang et al. 2006; Lopez-Sepulcre et al. 
201 1; Takahashi & Ho 2012), the overall picture suggests that 
there are two outflows in the region, as already proposed by 
Sanchez-Monge (2011): outflow A, most likely driven by a 
Class protostar; and outflow B. Thus, the CO (2-1) knots 
observed with the PdBI-A could be mainly tracing the cavity 
walls of outflow B, as most of the knots are seen in pairs along 
the direction of the HCO + high-velocity emission, and the 
strong CO (2-1) blueshifted knot 5" to the southwest of MM2 
(as well as the CO (2-1) redshifted emission to the northeast) 
could be associated with outflow A. MM2, which is a Class 0- 
like source, not detected up to 4.5 fim, and with emission from 
complex organic molecules (Palau et al. 201 1), is most likely 
the driving source of outflow A. As for the driving source of 
outflow B, this could be either an unresolved close companion 
of MM2 (at < 300 AU) or MM2-S itself (600 AU to the south 
of MM2), which falls within < 1" of the only 3.6 yum Spitzer 
source detected in the field. Therefore, in this region there is 
one strong source dominating the millimeter emission, which 
is driving at least one outflow. 

A5142: — The region has been proposed to harbor a small 
protocluster, because different millimeter sources (Zhang et 
al. 2007), at least three outflows (Hunter et al. 1995, 1999; 
Chen et al. 2005, Zhang et al. 2007) and multiple water masers 
(Hunter et al. 1995, 1999) have been reported. Our PdBI-A 
observations reveal a striking richness of millimeter sources 
(especially compared to the other regions), with two strong 
sources in the center and 5 additional sources all found within 
7" (12000 AU or 0.06 pc), as well as multiple CO (2-1) high- 
velocity knots (Fig. O. The strongest CO (2-1) redshifted 
knots are all located along the northwest direction following 
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RA (J2000) RA (J2000) 

Fig. 5. — Left: Same as Fig.f5]for A5142 (grey contours correspond to the integrated emission of NH3(1,1), Zhang et al. 2002). Contours of the PdBI CO (2-1) 
emission start at 55% and increase in steps of 10% of the maximum, 1 .85 Jy beam -1 km s _I , for the blueshifted emission, and start at 45% increasing in steps of 
10% of the maximum, 3.18 Jy beam -1 km s , for the redshifted emission. CO (2-1) blueshifted emission has been integrated from -19.7 to -8.9 km s , while 
CO (2-1) redshifted emission has been integrated from 0.7 to 12.7 km s _I (systemic velocity is -3.9 km s~'). Stars mark the positions of the hot molecular cores 
(Palau et al. 2011), shown on the right panel; plus signs indicate the positions of 2MASS sources and crosses indicate the positions of the 3.6 /im Spitzer/IRAC 
sources. Right: zoom on the central 4" of A5142. Thin black contours correspond to the 1.3 mm emission shown on the left panel. Thick green contours 
correspond to emission of (CH3)2CO tracing the hot molecular cores (Palau et al. 201 1), labeled MM1-HC (northern) and MM2-HC (southern). The synthesized 
beam of the PdBI CO (2-1) emission, of 0.56" X 0.45" with P.A.= 24.9°, is shown in the bottom-right corner. 



the well-known outflow C (e. g., Zhang et al. 2007), and the 
fainter redshifted knots seen to the southwest are well aligned 
with water masers (Goddi & Moscadelli 2006) and CO emis- 
sion from outflow B (Zhang et al. 2007). Additionally, there 
are hints of a possible new outflow (knots 'D'), which would 
consist of 3 blueshifted knots well aligned with MM2 and 
MM6, and extending roughly along the east-west direction. 
Hints of such an east-west outflow were also found in the 
blueshifted SO emission by Zhang et al. (2007, see figure 7 
of that paper). The driving sources of the outflows seem to be 
all located at the central 2" of the protocluster, where MM1 
and MM2 are found together with two additional hot cores, 
MM1-HC and MM2-HC, detected in different complex or- 
ganic molecules (Fig. fright, Palau et al. 201 1). As explained 
in Palau et al. (201 1), MM1-HC is the driving source of out- 
flow C, as its peak falls at the center of symmetry of a thermal 
radiojet elongated in the southeast-northwest direction (Zhang 
et al. 2007; Goddi et al. 2011), and it is also coincident with 
the center of symmetry of CH3OH masers tracing an infalling 
disk (Goddi et al. 2011). Concerning MM1, its millimeter 
continuum emission elongated in the east-west direction sug- 
gests that this source could be driving the north-south outflow 
A. Outflow B and the blueshifted knots 'D' could be related 
to MM2 and/or MM2-HC, as already proposed by Goddi & 
Moscadelli (2006) from water maser observations (note that 
MM2 is elongated in the north-south direction, thus perpen- 
dicular to knots 'D\ and MM2-HC is elongated perpendicular 
to outflow B, see Fig. |5}right). As for the 5 additional mil- 
limeter sources surrounding the central part of the protoclus- 
ter, MM3, MM6 and MM9 are found exactly along outflow 
axes, and it is not clear if they are tracing dust entrained by 
the outflows, as recently found by Maury et al. (2010, 2012), 
or are true protostellar objects. The fact that MM3 and MM6 
are associated with water maser emission at velocities about 



~ 5 km s away from the systemic velocity would rather 
suggest that they are fragments entrained by the passage of 
the outflow. The nature of MM7, close to a 2MASS star, and 
MM8, not exactly falling along any outflow axis, could be 
protostellar. In summary, this region presents a rich variety 
of millimeter sources, with two of them located at the center, 
dominating the millimeter emission and which seem to drive 
outflows, and the other fainter millimeter sources being either 
fragments entrained by the passage of the outflows, or true 
protostellar deeply embedded low-mass sources. 

4. ANALYSIS 

4.1. Regions in the literature used to complete the sample: 
general properties of massive dense cores 

In order to enlarge our sample, we searched the literature 
for massive dense cores observed with millimeter interferom- 
eters down to spatial resolutions of ~ 1000 AU and with mass 
sensitivities < 0.5 M @ . We required such a spatial resolution 
and sensitivity because this is the resolution and sensitivity of 
the optical/infrared studies reporting a correlation between the 
richness of the cluster and the mass of its most massive star 
(e. g., Hillenbrand 1995; Testi et al. 1999), and our final goal 
is to study such a correlation in the millimeter/sub millimeter 
ran ge. The final list, of 18 regions in total, is given in Ta- 
ble EE3, where we list the main observational parameters for 
each region. The rms noise of each image, as well as the 
angular resolution, field of view, and largest angular scale to 

14 In our literature search we found two regions, Orion-KL (Beuther et al. 
2004) and G5.89 (Hunter et al. 2008), which fulfill our selection criteria but 
were not included for different reasons. Orion-KL is a region undergoing a 
dynamical interaction between multiple stars (e. g., Zapata et al. 2009, 201 1; 
Moeckel & Goddi 2012), which is a very rare case in the literature; and the 
1 mm continuum emission of G5.89 is strongly contaminated by free-free 
emission (Hunter et al. 2008). 
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TABLE 3 

Massive dense cores studied with interferometers at 1.3mm with high sensitivity and down to a spatial resolution < 1000 AU 
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1.80 


2200 


55 
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a T in front of a source name stands for 'IRAS'. References: 1: Neri et al. (2007); 2: this work; 3: Carrasco-Gonzalez et al. (2012); 4: 
Fuente et al. (2005); 5: Fuente et al. (2007); 6: this work; 7: Zapata et al. (2005); 8: this work; 9: Leurini et al. (2007); 10: Cesaroni et 
al. (2005), Moscadelli et al. (2011); 11: this work; 12: Fernandez-Lopez et al. (2011); 13: Rodon et al. (2008); 14: Wang et al. (2012a), 
Jimenez-Serra et al. (2012); 15-18: Bontemps et al. (2010). 

b Bolometric temperature, defined as the temperature of a blackbody having the same mean frequency as the observed spectral energy 
distribution (e. g., Chen et al. 1995). 

c Rms of the millimeter images compiled for this work (see tablenote 'a' for references). M min is the mass sensitivity derived assuming 5 
times the rms noise, a dust temperature of 20 K, and a dust mass opacity coefficient at 1.3 mm from Ossenkopf & Henning (1994) as in 
Table [2] 

d The spatial resolution is calculated from the synthesized beam of the millimeter images compiled for this work (see tablenote 'a' for 
references), and the adopted distance for each core. LAS is the largest angular scale detectable by the interferometric observations, 
estimated following Palau et al. (2010) and converted to spatial scale using the adopted distances. 

e Mm,, is defined as the number of millimeter sources within ~ 0.1 pc (in diameter, see Section l4~2l >. N m m+m = Nmm + Mr, with N\ R 
being the number of infrared sources (2MASS + Spitzer/IRAC at 3.6 fan) without millimeter counterpart within ~ 0. 1 pc. In case of 4 
or 5o" detection but with independent claims in the literature of being a true source, we counted the source has 'half (0.5). Positional 
uncertainties are adopted as ~ 2" for Spitzer/IRAC, and ~ 1" for 2MASS and millimeter sources, and two sources were considered to 
be associated if their angular separation is smaller than the sum of the positional uncertainties. 

f First approach to (proto)stellar density estimated from N mm+ m divided by the volume of a sphere with the same projected area onto the 
sky as the polygon with minimum area (A po i y ) encompassing all the sources. For the case of iV mm+ iR < 2, we used the average spatial area 
measured for regions with A^m+iR > 2. The separation is an 'average' projected separation between sources in each region, estimated as 

-y^poly/Wmm+IR- 

g Observations of 3 mm continuum emission, and we assumed a spectral index of 3 to estimate the 1 mm flux densities. 



which the interferometer is sensitive have been converted to 
mass sensitivity and spatial scales to ease the comparison be- 
tween regions. The last 4 regions of Table [3] correspond to 
4 regions studied by Bontemps et al. (2010) which fulfill our 
selection criteria (Section 1). These 4 regions are separated 
from the others because they were observed with the PdBI in 
CD configuration, recovering spatial structures about a factor 
3-5 larger than the other regions of Table [3] (see column 9 in 
Table [3). A compilation of the millimeter high angular reso- 
lution images of all the regions (for the last four see figures 2 
and 3 of Bontemps et al. 2010) are presented in Figs. [TTI and 
[T2"lof the Appendix. 

In order to characterize in a uniform way the main proper- 



ties of the massive dense cores of Table [3j we have analyzed 
with the same methodology (using the flux densities reported 
by Motte et al. (2007) and di Francesco et al. (2008) for each 
core, and calculating all the parameters as described below) 
the following parameters for all the massive dense cores of 
Table [3] (results listed in Tables[3]and[4]i. 

Bolometric luminosity and bolometric temperature (Lj, i, ?boi) — 
To estimate Lboi and 7boiE3, we compiled the flux densities 
at different frequencies from 2MASS, Spitzer/IRAC+MIPS, 

15 is denned as the temperature of a blackbody having the same 
mean frequency as the observed spectral energy distribution (e. g., Chen et 
al. 1995). 
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TABLE 4 

Properties of massive dense cores studied with interferometers at 1.3 mm with high sensitivity and down to a spatial resolution < 1000 AU 
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200 


44 


31800 


2.25 


1.9 


4.3 


6 


1.3 


0.37 


16 


0.004 


11-I22134-IRS1 


250 


47 


65300 


0.67 


0.3 


0.1 


0.5 


3.4 


0.33 


34 


0.072 


12-HH80-81 


360 


60 


32900 


3.80 


3.1 


9.3 


9 


1.0 


0.67 


51 


0.041 


13-W3IRS5 


2500 


56 


59100 


8.03 


3.6 


2.2 


> 0.4 


0.9 






0.021 


14-AFGL2591 


1300 


140 


83300 


2.19 


0.7 


3.9 


0.9 


2.1 






0.013 


15-CygX-N53' 


68 


4 


19700 


1.98 


2.7 


< 2.2 


< 27 


1.1 








16-CygX-N12' 


70 


5 


16800 


2.79 


4.4 


< 1.3 


< 16 


0.8 








\l-CygX--H6y 


47 


10 


10100 


5.18 


14 


< 6.0 


< 60 


0.5 








18-CygX-N48 i 


160 


28 


22500 


3.54 


4.2 


< 1.3 


< 8 


0.9 









a T in front of a source name stands for 'IRAS'. 



b Mass of the dense core measured with the JCMT at 850 fim or the IRAM30m at 1.2 mm (Jenness et al. 1995; Motte et al. 2007; di 
Francesco et al. 2008), for a dust temperature of 20 K, and a dust (+gas) mass opacity coefficient of 0.0175 cm 2 g" 1 at 850 /im and 
0.00899 cm 2 g" 1 at 1.2 mm (Ossenkopf & Henning 1994). 

c Core size is the deconvolved FWHM. Surface density (S) and mean particle density (n) of the dense core are estimated assuming 
spherical symmetry and using X = M SD /(nR 2 ) (with R = FWHM/2) and n = 3 M SD /(4^ FWHM) (Motte et al. 2007). 

d ^max is the mass of the strongest millimeter condensation measured by an interferometer with a very extended configuration (see Table[3] 
for the references), assuming a dust temperature of 50 K and the dust opacity law of Ossenkopf & Henning (1994) as in Table|2] The 
uncertainty in the mass estimate due to the temperature assumption is a factor of 2. 

e The Core Formation Efficiency is calculated from the ratio of the sum of the masses of all the fragments measured with an interferometer 
in very extended configuration (Table O and the single-dish mass derived from the peak intensity, making the same assumptions as 
tablenote 'b\ Lower limits correspond to single-dish masses possibly contaminated by other nearby cores. 

f M jLs is the thermal Jeans mass for a temperature of 20 K (M* ans = 0.6285 [r/10K] 3/2 [;7 H2 /10 5 cirr 3 F 1/2 M a ). M\f ms is the Jeans mass 
estimated using the total velocity dispersion (thermal+non-thermal, Pillai et al. 201 1). 

g CTnon-th is the non-thermal velocity dispersion, calculated by fitting the hyperfine structure of NH 3 (1,1) emission of a quiescent core in the 
surroundings (< 2') of our studied massive dense core (see references in Section l4~71 (. and assuming a thermal dispersion corresponding 
to 20 K. 

h /3 rot : ratio of rotational-to-gravitational energy, as defined in Chen et al. (2012, see Sect. 14. lb . derived from dense gas tracers tracing the 
massive dense core at scales of 5"-10". 

1 Regions from Bontemps et al. (2010) observed with the PdBI in CD configuration. This implies that some of the parameters are upper 
limits when compared to the rest of the sample (observed with more extended configurations). The parameters are essentially the same 
as reported in Bontemps et al. (2010), using the opacity law adopted in this work and the new distance of Cygnus X, 1.4 kpc (Rygl et al. 
2012). 



WISE, MSX, IRAS, and JCMT (SCUBA, di Francesco et 
al. 2008) catalogs. This allowed us to build the spec- 
tral energy distribution (SED) of the (presumably) most lu- 
minous source embedded in the massive dense core. For 
some regions where Spitzer data were strongly saturated, 
we used FIR/submillimeter fluxes from the literature (But- 
ner et al. 1990; Campbell et al. 1995). When possi- 
ble, we compared our compiled photometry (especially the 
Spitzer/IRAC+MIPS and WISE fluxes) with fluxes reported 
in the literature at similar wavelengths (e. g., Lada et al. 1984; 
Skinner et al. 2009; Crimier et al. 2010), and found them to 
be consistent. To build the SED in the infrared we used the 



flux densities of the 2MASS and Spitzer/IRAC source clos- 
est to the strongest millimeter source detected with an inter- 
ferometer, which fell within < 5" of the Spitzer/MIPS and 
JCMT/SCUBA peaks of emission, for all the cases except for 
OMC-1S. [3- Thus, the compiled photometry was derived 
from data taken with beams < 20", similar to our studied 

16 For the case of OMC-1S, the infrared source chosen to build the SED is 
source 136-355 (Zapata et al. 2005), which is about ~ 12" to the north of the 
strongest mm source. Given the number and intensity of millimeter sources 
embedded in the dense core traced by the JCMT at 450 /im (di Francesco et 
al. 2008), we estimated that about 1/3 of the total submm flux density at 450 
and 850 /im comes from source 136-355, and used these estimates to build 
the SED of this source. 
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10000 100000 10 6 

n H2 (cm" 3 ) 

Fig. 6. — CFE vs mean density of the core, «h, , after Bontemps et al. (2010). 
Green triangles: CygnusX cores (Bontemps et al. 2010); grey triangles: p- 
Oph cores (Motte et al. 1998); blue dots: cores presented in this work. For 
Bontemps' and Motte's cores, the CFE is an upper limit (compared to the 
other cores) because the values are measured with a compact configuration 
of an interferometer. 

field of view. The obvious exception are the IRAS fluxes, but 
IRAS, Spitzer/MIPS and WISE were found to be consistent 
within 30^-0%. This indicates that the IRAS fluxes have their 
main contribution from the strongest millimeter source of the 
field and that the effects of multiple sources within the IRAS 
beam are not dominant in our sample. Overall, we estimated 
an error of about 30^40% in the derived bolometric luminos- 
ity of all the regions, except for OMC-1S, for which the error 
should be of a factor of 2 because its SED is the only SED not 
properly sampled between 3.6 /mi and 450 //m. Therefore, 
we are confident that the Lboi given in Table [3] is a reason- 
able approach to the luminosity (and hence to the mass) of the 
most massive star embedded in each massive core. Concern- 
ing 7b i, our derived values span a range of one order of mag- 
nitude, from 30 to ~ 200 K, with 10 cores having T\,oi< 60 K. 

Mass, mean surface density, and mean density of the core ( Msd. S, 
n H J — The mass of each dense core was estimated using the 
JCMT/SCUBA flux densities at 850 fim (di Francesco et al. 
2008, Jenness et al. 1995), or the 1 .2 mm flux densities for the 
Cygnus X regions (Motte et al. 2007), adopting the new dis- 
tance to CygnusX of 1.4 kpc (Rygl et al. 2012). To estimate 
the mass we assumed a dust temperature of 20 K, and a dust 
(+gas) mass opacity coefficient of 0.0175 cm 2 g -1 at 850 fim 
and 0.00899 cm 2 g 1 at 1.2 mm (Ossenkopf & Henning 1994). 
The masses range from 17 to 2500 M e , and increase with Lboi 
(in Table |U the regions are ordered with increasing luminos- 
ity), as already seen in previous works (e. g., Sridharan et al. 
2002). By using the FWHM of the 850 fim core (measured in 
the images of di Francesco et al. 2008, and reported in Motte 
et al. 2007), and assuming spherical symmetry, we derived 
the surface density and mean densitvPl of each core. These 
parameters range from 0.7 to 10 gem 2 , and from 0.3 to 30 
times 10 5 cm 4 . 

Ratio of bolometric luminosity to mass of the core (L^qi/Msd) — 
Lboi/MsD was estimated as an alternate measure (together with 
7boi) of the evolutionary stage of the dense core (e. g., Sridha- 
ran et al. 2002; Lopez-Sepulcre et al. 2011). While L\, i/Msd 
is more sensitive to Msd and thus more representative of the 

17 We estimated the surface density following Z = M/(nR 2 ) (with R = 
FWHM/2) and the mean density following n = 3 M/(4tt FWHM), as in 
Motte et al. (2007) and Bontemps et al. (2010). 




1000 10000 100000 



Lboi(Ls) 

Fig. 7. — M milx vs L[, i for the (sub)sample of 14 sources. The figure shows 
a possible weak trend between M milx (the mass of the most massive fragment 
as measured with an interferometer in extended configuration) and L(, i- In 
this plot 122134 is not included because the two most massive sources of this 
region (t he U CHII region and IRS1), are almost completely deprived of dust 
(Section 13. 2t . and cannot be compared to figure 2 of B onne ll et al. (2004), 
calculated at the time of protostar formation (see Section |5~j"l . 

dense core as a whole, 7b i is more sensitive to the proper- 
ties of the infrared source chosen to build the SED, as Lboi 
assesses the mean frequency of the SED weighted with the 
flux density at each frequency. In the following, we will use 
Lboi/MsD as a proxy to the core evolutionary stage. 

Core Formation Efficiency (CFE) — We calculated the CFE as 
the ratio between the sum of the masses of all the fragments 
measured with an interferometer in very extended configu- 
ration (Table [3) and the single-dish mass estimated from the 
peak intensity. We used the single-dish peak intensity because 
the interferometer field of view, ~ 20", is similar to the single- 
dish beam (~ 23" at 850 fim, di Francesco et al. 2008). The 
values of CFE range from 0.4 to 60%. Interestingly, we find a 
trend (Pig. [6]) between CFE and the mean density of the core, 
which fits within the trend previously found by Motte et al. 
(1998) and Bontemps et al. (2010; see grey and green dots in 
Fig. [6]), and which is suggestive of gravity playing a dominant 
role in the densest cores, favoring the formation of compact 
structures. 

Mass of the strongest millimeter source ( M raax ) — M mdx was cal- 
culated from the flux density of the strongest millimeter 
source measured by an interferometer with a very extended 
configuration (see Table [3] for the references; the CygnusX 
sources have been excluded as they were observed with the 
PdBI in CD configuration), assuming a dust temperature of 
50 K and the dust opacity law of Ossenkopf & Henning (1994) 
as in Table 12 Some of these sources are associated with hot 
molecular cores for which temperatures of ~ 100 K have been 
estimated. However, this information is not available for all 
the sources, and for the temperature range of 30-100 K the 
masses change by a factor of 2 with respect to the value de- 
rived assuming 50 K, which we consider a reasonable first 
approach. The plot of M max vs Lboi (Fig. shows a we ak 
correlation, which will be further discussed in Section lBTTI 

Non-thermal velocity dispersion (o- non _ th ) — To have a first idea 
of the initial non-thermal velocity dispersion in the core, we 
searched for interferometric data of NH3(1,1) reported in the 
literature and found data for 8 regions (Zhou et al. 1990; 
Tieftrunk et al. 1998; Zhang et al. 1998; Zhang et al. 2002; 
Gomez et al. 2003; Busquet 2010; Keto & Zhang 2010; 
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Fig. 8. — N mm vs a number of properties etimated for the 18 massive dense cores and compiled in Tables F5] and f4] (see Section PTTV The numbers on the top 
of each dot indicate the source-ID listed in TablesfJ]and|4] Blue dots correspond to sources labeled from 1 to 14 in the tables, and the green dots correspond to 
sources 15 to 18, which are the Cygnus X sources studied in Bontemps et al. (2010). In panel H, the vertical dotted line marks the sound speed (c s , corresponding 
to a temperature of 20 K). 



Sanchez-Monge 201 1; Sanchez-Monge et al. in prep.; all ob- 
served with the Very Large Array with a typical angular res- 
olution of ~ 5" and field of view ~ 2')- Among the different 
NH3 cores surrounding the massive dense core of our sample, 
we chose the most quiescent core, which should be barely 
affected by turbulence injected by outflows or by systematic 
motions, and should be representative of the initial turbulence 
of the massive dense core. With the measured linewidth in this 
quiescent core, we estimated the non-thermal velocity disper- 
sion, cr non _ t h (assuming a thermal component corresponding 



to 20 K), and found values from 0.3 to 0.7 km s _1 0, all su- 
personic (sound speed corresponding to 20 K is 0.27 km s ). 
The narrow range of values estimated for <x n on-th, around 
0.5 km s _1 , suggests that our approach is reasonably avoid- 
ing outflows/infall/rotation motions. 

18 For the Cygnus X sources, Bontemps et al. (2010) estimate the non- 
thermal velocity dispersion from single-dish measurements at the position 
of each massive core, with an angular resolution of ~ 20", thus possibly 
including turbulence injected by outflows/systematic motions, and for this 
reason we do not include these values here. 
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Jeans mass and Jeans number (Mj Mns , Nj tllm ) — We estimated 
the 'mean' Jeans mass using the mean density of the cores 
listed in Table|4j and assuming a gas temperature (as in gravo- 
turbulent fragmentation scenarios, e.g., Schmeja & Klessen 
2004; Dobbs et al. 2005), and found M Jeans -0.5-3 M Q (Ta- 
ble 0J. This is much smaller than the mass of the core, and 
larger than the mass sensitivity achieved for regions of our 
sample (Table [3}, implying that we should have detected be- 
tween 20 and 200 fragments if a pure thermal Jeans frag- 
mentation was at play. We also calculated Mj eans using the 
total velocity dispersion (thermal+non-thermal, also called 
'turbulent Jeans' fragmentation, Pillai et al. 2011), yielding 
Mjeans SS M^d, which gives probably an idea of the mass 
of the most massive fragments expected. Thus, according 
to Jeans fragmentation, one would expect a large number of 
low-mass fragments and a few massive fragments, while we 
observe only a few low-mass fragments. 

Rotational-to-gravitational energy (j3 mt ) — Simulations show 
that rotating cores fragment if their rotational energy is large 
enough compared to the gravitational energy (e. g., Boss 
1999; Forgan & Rice 2012). The threshold derived from the 
simulations is /3 rot > 0.01, and recent observations of a sample 
of low-mass rotating systems seem to be consistent with this 
lower limit (Chen et al. 2012). In our sample, [5 ml was esti- 
mated following equation (1) of Chen et al. (2012), and using 
interferometric observations of NH3(1,1) (for 8 regions, see 
above), NH 3 (2,2), CS(2-1), and S0 2 (Torrefies et al. 1989; 
Beltran et al. 2004; Rodon et al. 2008), tracing the core at 
scales of 5-10". We caution that the velocity gradients seen 
in the dense cores, used to estimate f3 rot , could be affected in 
some cases by multiple velocity components and/or perturba- 
tion by outflows. Despite of this, the estimated values of /3 rot 
range from 0.003 to 0.072 (Table |4j, similar to the range of 
values reported in the literature (e.g., Goodman et al. 1993; 
Chen et al. 2012). 

4.2. Fragmentation level in massive dense cores 

The total sample of 18 massive dense cores compiled in Ta- 
ble [3] spans a range of bolometric luminosities from ~ 300 
up to 2 x 10 5 Lq, or almost 3 orders of magnitude, and we 
aim at relating the dense core properties, such as Lboi, with 
their fragmentation level. To achieve this, we did the follow- 
ing. For each high angular resolution millimeter image (see 
Table [3] for references) we estimated the rms noise using the 
same procedure in Mapping of the Gil das s oftware and identi- 
fied those sources above 6<x (in Section lXTI we showed that the 
6<x cutoff yields reliable identifications; the 6cr level is marked 
in red in the cleaned images of Figs. [TTI and [T2l. We defined 
A^mm as the number of millimeter sources inside a region of 
0. 1 pc of diameter (marked in the figures of the Appendix as 
a black circle). In case of 4 or 5 cr detection but with indepen- 
dent claims in the literature of being a true source, we counted 
the source as 'half (0.5). The median number of N mm in the 
sample is 3.25. 

Among the 18 studied cores in Table[3] five (~ 28%) show 
a very low degree of fragmentation, with only one source 
(or one dominant source with a much fainter companion) 
within the 0.1 pc region, while 9 (~ 50%) cores split up into 
> 4 sources. The sources with a low degree of fragmen- 
tation are prototypical deeply embedded intermediate/high- 
mass sources (e.g., Cesaroni et al. 2005; van der Tak et al. 
2006; Fuente et al. 2005; Sanchez-Monge et al. 2010; Bon- 
temps et al. 2010), as their low 7boi and Lboi/A^SD indicate (see 



Tables [3] and Therefore, the observed low fragmentation 
is not an evolutionary effect of clusters which have already 
cleared out their natal material. 

In addition, in order to have a first estimate of the infrared 
sources possibly associated with each massive dense core, 
we searched the 2MASS Point Source Catalogue (Skrutskie 
et al. 2006) and downloaded the images from Spitzer/IRAC 
databas43 to extract the positions of the 3.6 jum sources 
within the 0.1 pc field of view for each core. The final posi- 
tions of the infrared (2MASS and 3.6 fim) sources are marked 
in Figs.[TT]and[T2] We counted the number of infrared sources 
(2MASS + Spitzer/IRAC at 3.6 /urn) within 0.1 pc (see Ta- 
ble O and found that the total number of millimeter and in- 
frared sources, NiR +mm , correlates with N mm , i. e., the regions 
with a small number of millimeter sources also seem to have a 
small number of infrared sources, and vice versa (see Table[3] 
and Fig.[8]-d). Similarly, the ratio Nig./N mm is in general close 
to 1 (TableO. We are aware that there are several factors con- 
tributing to the uncertainty of /Vir. On the one hand, the mass 
sensitivity of 2MASS is not as good as the sensitivity of the 
PdBI-A, especially if there is important visual extinction in 
the region (as is the case for our regions). We estimated that, 
following the method of Testi et al. (1998), for an (optimistic) 
visual extinction of Ay ~ 5 mag, the mass sensitivity of 
2MASS in A^-band (which has a limiting magnitude of 14.3, 
Skrutskie et al. 2006) is 0.6-1.4 M Q , for an age of 0.5-1 Myr. 
On the other hand, we did not do any color-color analysis to 
discard possible foreground/background contaminants. As for 
the Spitzer/IRAC 3.6 fim data, for a region at 2 kpc, we are 
sensitive to ~ 0.1 M (using a limiting magnitude of ~ 16.5, 
Qiu et al. 2008). The only caveat with Spitzer/IRAC images 
is that the two brightest sources of the sample (W3IRS5 and 
AFGL2591) are affected by saturation effects, and thus /Vir 
is a lower limit in these two cases. In spite of these short- 
comings, it is striking that the ratio N\R/N mm is — 1 for 12 
out of 18 sources and, although the exact number of infrared 
sources embedded in the massive dense cores should be stud- 
ied more accurately, the trend shown in Fig. [8}d is probably 
real. Overall, the total number of millimeter sources found 
in all the regions is 57, of which only 20 have infrared coun- 
terparts. This indicates that only ~ 30% of the detected mil- 
limeter sources are also detected in the infrared, consistent 
with previous studies by Qiu et al. (2008), and indicating that 
in general the millimeter sources in our sample of massive 
dense cores are deeply embedded and in very early evolution- 
ary stages. 

As a last calculation, we assessed the protostellar density by 
dividing the total number of millimeter+infrared sources (Ta- 
bleO by the volume corresponding to a sphere with the same 
projected area onto the sky as the polygon with minimum area 
encompassing all the sources. The result is shown in Table[3] 
togethter with the mean separation between sources, typically 
of ~ 3000 AU. These numbers should be regarded as prelim- 
inary and as lower limits for W3IRS5 and AFGL2591 due to 
limitations of the infrared datasets (see above). The average 
stellar density is around ~ 10 5 stars pc -3 , similar to the stellar 
density measured in the Orion Nebula Cluster (e.g., Hillen- 
brand & Hartmann 1998). 

5. DISCUSSION 

19 Data downloaded fr om the 2MASS and Spitzer Her- 
itage Archive webs, http://www.ipac.caltech.edu/2mass/ and 
http://sha.ipac.caltech.edu/applications/Spitzer/SHA. 
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5.1. Implications for intermediate/high-mass star formation 

In Section |4~Tl we found a weak trend between the mass of 
the strongest millimeter source and the bolometric luminos- 
ity estimated for each core (Fig. [7). In the original competi- 
tive accretion scenario, Bonnell et al. (2004) predicted that the 
mass of the small-scale (1000-2000 AU) condensation asso- 
ciated with the nascent intermediate/high-mass star does not 
correlate with the final stellar mass at the time of protostar for- 
mation (see figure 2 of Bonnell et al. 2004). The weak trend 
shown in Fig.[7]does not support a pure competitive accretion 
scenario. On the other hand, Fig. [6] shows that for most of 
the cores only a small fraction of mass is found in compact 
fragments, which means that most of the mass is in an ex- 
tended envelope outside the fragments. This finding counters 
the monolithic collapse scenario (McKee & Tan 2002, 2003). 
Thus, the results presented here (Figs.[6]and[7]) show that nei- 
ther monolithic collapse nor competitive accretion accurately 
prescribe the physics in these regions. Zhang & Wang (201 1) 
propose an alternative picture that condensations/fragments 
continue to accrete mass from the general environment (called 
'core' in this paper), whereas the embedded protostar accretes 
mass from the condensation/fragment. This 'simultaneous ac- 
cretion picture' explains both Figs. [6] and |7j 

5.2. What is controlling fragmentation in massive dense 
cores? 

In the previous section we presented a sample of 18 massive 
dense cores observed in the millimeter range down to a spa- 
tial resolution of ~ 1000 AU and mass sensitivities < 0.5 M . 
Among the 18 cores, about ~ 30% show no fragmentation, 
while 50% harbor > 4 fragments. Given the different frag- 
mentation levels in our sample, we consider in the following 
possible ingredients which could play a role in the fragmen- 
tation of a core. 

Optical and infrared studies suggest that Lt, i is related to 
the number of sources in a cluster (e.g., Hillenbrand 1995; 
Testi et al. 1998). We searched for any correlation between 
N mm and Lt, i, and the result is plotted in Fig. [8}a. We find 
no signs of any possible correlation: the highest values of 
N m m are reached in regions of low/intermediate luminosity, 
such as OMC-1S and A5142, and regions with similar lu- 
minosity show very different properties (compare for exam- 
ple NGC7129-FIRS2, Fuente et al. 2005, to IC1396-N, Neri 
et a l. 20 07). Since Lboi correlates well with Msd (see Sec- 
tion I4.lt . this could indicate that the mass of the core is not 
a crucial ingredient in determining its fragmentation (see also 
Fig.|8}c). This could also indicate that the correlation seen in 
the optical/infrared is probably not established at early stages 
of cluster formation, and that the final number of members in 
a cluster is set after several star formation episodes, as sug- 
gested by both theoretical (e.g., Myers 2011) and observa- 
tional (e.g., Williams et al. 2009; Hennemann et al. 2010; 
Bik et al. 2012) works. A deeper study including a larger 
sample and larger field of view should be carried out to draw 
firm conclusions on this. We note however that cluster evo- 
lution alone is not likely enough to explain the different frag- 
mentation, since for cores with a similar value of the evolu- 
tionary indicators (7b i, L\, \/M^) we already find very dif- 
ferent fragmentation levels (see Fig. |8}b). This was in fact 
also suggested by the trend found between N mm and A^mm+iR 
(Fig.|8}d; or Ni R /N mm ~ 1, Table indicating that the differ- 
ent episodes of star formation yield comparable numbers of 
fragments, with some cores highly fragmenting from the be- 



ginning and keeping this property constant with time, while 
others remaining unfragmented during all core evolution. 

The initial density profile of a core could be one of the 
agents determining its fragmentation level, as simulations 
show that highly concentrated cores (e.g., «h, °c r~ 2 ) lead 
to the formation of one single massive star, while flat density 
profiles («h 2 k r or flatter) produce an important number 
of low-mass stars (e. g., Myhill & Kaula 1992; Burkert et al. 
1997; Girichidis et al. 2011). In our sample, two of the 5 
cores showing no fragmentation present concentrated density 
profiles (Beuther et al. 2002a; Mueller et al. 2002; Crimier et 
al. 2010), supporting this view. A uniform study of the density 
profiles of the cores presented in this work would be required 
to accurately assess its importance on the fragmentation level. 

Another ingredient which could be related to fragmentation 
is the rotational-to -gra vitational energy, yS rot . The calculated 
values of j3 mt (Sec. 14. II Table @]i present no relation with N mm 
(Fig. [8}f), although on average, the low-fragmented regions 
have /Jrot smaller than the highly fragmented regions. The 
highly fragmented regions present a wide range of f> rot values 
spannig at least one order of magnitude. This could indicate 
that y8 rot is not an accurate measure of the rotation of the cores 
(i. e., the velocity gradients measured may not always trace 
rotation), or that a wide range of angular momenta can yield 
similar fragmentation. 

The turbulence associated with a core could be also deter- 
mining its fragmentation. For example, 'turbulent fragmenta- 
tion' models show that turbulence develops shocks which in 
turn produce substructure in the core (e. g., Padoan & Nord- 
lund 2002; Klessen et al. 2005), and one could expect a 
higher fragmentation for higher turbulence levels (Jappsen & 
Klessen 2004; Peretto et al. 2007). However, N mm does not 
seem to depend either on the non-thermal (tur bulen t) velocity 
dispersion of the core s of o ur sample (Section l4"Tl Fig.[8}h). 

Finally, in Section 14.11 we showed that there is a deficit 
in the observed number of low-mass fragments as compared 
to expectations from thermal Jeans fragmentation, and this 
is true even for cores with rather flat density profiles (e. g., 
122134, IRAS05358+3543NE, Beuther et al. 2002a). This 
is very similar to the conclusions drawn by Bontemps et al. 
(2010) and van Kempen et al. (2012) for their (samples of) 
sources. 

All this suggests that still other forms of support additional 
to thermal and turbulence may play a role. Interferometric 
observations of polarized submillimeter emission have shown 
that the magnetic field seems to play an important role in the 
dynamics of some massive dense cores (e. g., Lai et al. 2003; 
Girart et al. 2009; Tang et al. 2009). In the following section 
we discuss the influence of magnetic field and radiative feed- 
back on the fragmentation process of a massive dense core. 

5.3. Role of combined magnetic field, radiative feedback and 
turbulence in the fragmentation process 

A number of studies comparing observations to simula- 
tions conclude that magnetic fields could be a key ingre- 
dient in the fragmentation process. For example, for the 
low-mass case, Maury et al. (2010) find that magneto- 
hydrodynamical models agree much better with their obser- 
vations. And for the intermediate/high-mass case, Peretto 
et al. (2007) find difficulties matching the observed masses 
and number of fragments with the results of hydrodynami- 
cal simulations, suggesting that an extra support such as pro- 
tostellar feedback or magnetic fields is at play. In order to 
assess the role of magnetic fields and radiative feedback in 
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mu=2 




TABLE 5 

Comparison of Commercon et al. (201 1a) simulations to PdBI observations 
of this work within a field of view of ~ 0.07 pc 
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-5000 5000 
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Fig. 9. — Direct output of Commercon et al. (201 la) simulations (column 
density maps) for the highly magnetized core model (top) and the low mag- 
netized core model (bottom). 



the fragmentation of a massive core we considered the radi- 
ation magneto-hydrodynamic simulations of Commercon et 
al. (201 la). These simulations include the effects of magnetic 
fields, turbulence, and take into account the radiative feedback 
from the accretion shock (i.e. shock on the first Larson core) 
and thermal emission (no radiative feedback from the accret- 
ing/forming protostars using pre-main sequence stellar evolu- 
tion subgrid models), down to a spatial resolution of 32 AU 
(low resolution runs, see fig. 2 in Commercon et al. 2011a). 
We note that the initial conditions adopted by Commercon 
et al. (2011a) correspond to a relatively fiat density pro- 
file, which should favor fragmentation (e. g., Girichidis et al. 
201 1). Commercon et al. (201 la) find that magnetic fields and 
radiation can combine to inhibit fragmentation in a massive 
core, and that this is particularly important if the initial mag- 
netic intensity is strong. In their simulations, the initial mag- 
netic intensity is set by the parameter /u = (M '/ <t>) / \M '/ <t>) cr i t 
(mass-to-flux over critical mass-to-flux ratio, Mouschovias 
& Spitzer 1976; Hennebelle et al. 2011), and the two ex- 
treme cases explored are = 2 (highly magnetized core) and 
fi = 130 (weakly magnetized core). By using the output of 
the adaptive mesh refinement code RAMSES (Teyssier 2002; 
Fromang et al. 2006; Commercon et al. 2011b), we compute 
column density maps (in the three different xy, xz, and yz pro- 
jections) for each value of fi from the low resolution runs of 
Commercon et al. (201 la). 



Parameter 



Simulations 



Observations 



MU2 - highly magnetized core 



M 



(M ) a 



5000 AU 
Mjooau (M Q ) a 
CFE b 

number of fragments 



4.7 
3.8 
81% 
1 



3.2 
1.9 

59% 
1 



MU130 - weakly magnetized core 



M t0t (Mr-,Y 

5000AU V' H o/ 


16 


47 


M 5 00AU (M Q f 


4.0,3.9,1.7 


8.9,6.4,0.94, 




0.77,0.47 


0.85,0.81,0.77,0.72 


CFE b 


68% 


41% 


number of fragments 


5 


7 


fragment separation (AU) 


2000-3500 


2000-5000 



a Masses estimated assuming a dust temperature of 30 K for 
all the fragments, and the same opacity as in Table[2] M 50 oau 
refers to the mass of each fragment as observed with the PdBI 
in A configuration. 

b Core Formation Efficiency, defined as the ratio between 
the total mass measured with the PdBI in A configuration 
(~ 500 AU scale) to the mass measured with an interferome- 
ter in a compact configuration (~ 5000 AU scales). Note that 
this definition is different from the one used in Section |4~T1 
Table!] and Figs.|6]and[8] 



In order to compare the output of the simulations of Com- 
mercon et al. (2011a) with our observational interferometric 
data we selected two regions with low and high fragmenta- 
tion as seen with the PdBI-A. For the low fragmentation case 
we selected 122198 to be compared to the /i = 2 (highly 
magnetized) simulated maps, while for the high fragmenta- 
tion case we selected A5142 to be compared to the fi = 130 
(weakly magnetized) simulated maps. The column density 
simulated maps used for the comparison are shown in Fig. [9] 
which reveal a striking intrinsic difference between the highly 
magnetized (jj = 2) and the weakly magnetized (fi = 130) 
case. By using the output of the adaptive mesh refinement 
(AMR) code RAMSES (Teyssier 2002; Fromang et al. 2006; 
Commercon et al. 201 lb), we compute dust emission maps 
using the RADMC-3D code0 (in the three different xy, xz, and 
yz projections) for each value of yu from the low resolution 
runs of Commercon et al. (2011a). For the conversion to 
Jypx 1 , the 3D AMR structure of RAMSES is directly loaded 
into RADMC-3D and we assume that the gas temperature com- 
puted in the radiation-magnetohydrodynamics calculations of 
Commmercon et al. (2011a) is correct and equals the dust 
temperature. We use the opa city law of Ossenkopf & Henning 
(1994) at 1.3 mm as in Sec. 13. II We assume the correspond- 
ing distances for 122198 and A5142, and put the coordinates 
of 122198 to the yt — 2 map, and the coordinates of A5142 to 
the yu = 130 map. Finally, we convolved the resulting sim- 
ulated maps in Jypx 1 by the Mu-coverage of PdBI-A used 
in this worliPl inverted the visibilities and cleaned the dirty 

20 http://www.ita.uni-heidelberg.de/~dullemond/software/radmc-3d/ 

21 For the convolution of the simulated maps by the w«-coverage of PdBI- 
A, we used the GILDAS task 'uv_model' , which generates a set of visibili- 
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Fig. 10. — Top: Observations of 122198 (left) and A5142 (right) with the PdBI at 0.4" angular resolution (black contours) and the PdBI/SMA at l"-2" (grey 
contours; left: 122198 with the PdBI at 3 mm, Sanchez-Monge 201 1; right: A5142 with the SMA at 1.3 mm, Zhang et al. 2007). The field of view corresponds 
to the same spatial scale for both regions (0.07 pc). Bottom: Simulations convolved (see footnote 21) with the Mo-coverage of the PdBI-A (black contours) 
and PdBI-CD (grey contours). Bottom-left panel are the simulations for the magnetized core case (fi = 2). Bottom-right panel are the simulations for the low 
magnetization case (ji = 130). Left panels contours: Black: -4, 4, 8, 16, 32, 64, 128, and 256 times 2(1) mjy beam -1 for observations (simulations); grey 
contours: -4, 4, 8, 16, 32, 64, 128, and 256 times 0.6(1) mJy beam -1 for observations at 3 mm (simulations at 1 mm). Right panels contours: Black: —4, 4, 8, 
16, 32, and 64 times 2.8(1.4) mjy beam - for observations (simulations); grey contours: -2, —1, 1, 2, 4, 8, 16, and 32 times 12(1.2) mjy beam -1 for observations 
(simulations). 



simulated maps as done with our real observations, allowing 
us to compare the fragments, at ~ 1000 AU scales, that the 
PdBI-A would have detected if our regions were similar to 
the simulated ones. The results are shown in Fig. [10] (black 
contours; left panels of the figure for 122198, and right panels 
for A5142). In order to compare also the morphology of the 
observed massive dense core with the simulated one at scales 
~ 5000 AU, we convolved the simulated Jypx -1 images by 
the the wy-coverage of PdBI in a typical CD configuration of 
the PdBI (grey contours of Fig.fTOrbottom). 

As shown in Fi g . [TOl-left-bottom. for the highly magnetized 
core Qi = 2) the simulations convolved with the PdBI-A uv- 
coverage reveal one single fragment at the center, which is 

ties at the im-baselines of the PdBI-A observations (by calculating the Fourier 
transform of the simulated image), allowing us to mimick not only the syn- 
thesized beam of our observations, but also the effects of filtering large-scale 
structures by the interferometer. 



dominating the millimeter emission and which is slightly re- 
solved. This is very similar to what we observe in 122198 
(Fig. [TOlTeft-top'). where one single slightly elongated frag- 
ment dominates the millimeter emission down to a mass sen- 
sitivity of ~ 0.1 M Q , Concerning the weakly magnetized core 
(fi = 130), Fig. [TOl-right-bottom shows that the core splits up 
into five fragments, with 3 of them dominating the millimeter 
emission. This is again similar to what we observe in A5142 
(Fig. ITUl-right-top). where the large-scale envelope (Zhang et 
al. 2007) splits up into 6-7 fragments for a mass sensitivity 
of ~ 1 M , with two massive fragments in the center. The 
qualitative agreement between our PdBI observations and the 
simulations of Commercon et al. (201 la) is remarkably good. 

In Table [5] we present a preliminary quantitative compari- 
son of our observations and the simulations of Commercon 
et al. (2011a). The table compares the mass of the core at 
~ 5000 AU scales to the mass predicted by the simulations, 
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as well as the number of fragments, mass of the fragments, 
CFEpl, and average separation between fragments. Again, the 
comparison between Commercon et al. (2011a) simulations 
and our observations are consistent. For the // = 130 case the 
simulated mass of the dense core at scales of ~ 5000 AU is 
lower by around a factor of 3 with respect to the observations, 
but this could be the consequence of the «u-coverage used 
to reproduce the large-scale core (grey contours in Fig. flQb . 
which was not exactly the same as the wy-coverage of the 
A5142 large-scale observations (taken with the Submil- 
limeter Array, Zhang et al. 2007). In addition, the fact that 
the calculations are stopped at an arbitrary time after the 
formation of the first hydrostatic core is also affecting the 
mass of the simulated core (Commercon et al. 2011a). 

Therefore, the different fragmentation level seen in our 
sample of massive dense cores could be a consequence of dif- 
ferent initial magnetic field strength in the cores. For com- 
parison, the magnetic field strength assumed in the ji — 2 
case is around 0.3-1 mG at the scales of the massive dense 
core (~ 3000 AU), which is comparable to the magnetic field 
strength derived from observations of star-forming regions at 
similar scales, of about 0.5 mG (e. g., Girart et al. 2006; Fal- 
garone et al. 2008; Goncalves et al. 2008; Frau et al. 2011). 
This indicates that the magnetic field strength required to in- 
hibit significantly the fragmentation in a dense core for the 
simulations of Commercon et al. (201 la) can easily be found 
in star-forming regions, and drives us to carry out further ob- 
servational studies to study the relation of fragmentation and 
magnetic field strength in our sample. 

6. CONCLUSIONS 

We conducted 1.3 mm continuum and CO (2-1) PdBI ob- 
servations, using its most extended configuration, towards 
four massive dense cores (I22172N, 122134, 122198, A5142), 
allowing us to study their fragmentation for mass sensitivi- 
ties ~ 0.3 Mq and spatial resolutions ~ 1000 AU, compa- 
rable to studies of clusters in the optical/infrared. We com- 
plemented our high angular resolution observations with ob- 
servations from the literature of 14 additional massive dense 
cores observed at similar mass sensitivities and spatial reso- 
lutions. For the 18 sources in the sample we estimated a set 
of parameters in a uniform way and assessed their role in the 
fragmentation process. Our main conclusions can be summa- 
rized as follows: 

- The four regions presented in this work show a variety 
of fragmentation levels. Towards I22172N we find one 
strong and elongated millimeter source, driving a small 
outflow, and surrounded by 3 fainter millimeter sources. 
The mass of the strongest millimeter source is estimated 
to be around 2-10 M , and its size is of 2300x 1000 AU. 
122134 reveals four low-mass < 1 M & fragments dis- 
tributed surrounding the UCHII, for which we did not 
detect millimeter continuum emission. 122198 is the 
core with smallest degree of fragmentation, with one 
single source dominating the millimeter emission and 
no additional fragments down to masses of ~ 0.1 M Q . 
The strong millimeter source, with a mass of about 
~ 2 M Q and a size of 500 x 300 AU, is driving at least 
one outflow. Finally, for the case of A5142 we find the 

22 The Core Formation Efficiency (CFE) is measured here as the total mass 
in fragments (at ~ 500 AU scales) versus the mass of the core at ~ 5000 AU. 



highest degree of fragmentation. The millimeter emis- 
sion consists of two strong elongated sources, both driv- 
ing outflows, and with masses of about 5 M Q , which 
are surrounded by 5 additional unresolved fragments of 
about ~ 2 M Q . 

- For the total sample of 18 massive dense cores, we 
estimate luminosities in the range 300 to 2 x 10 5 L Q , 
single-dish masses of ~ 20-2000 M , mean densities of 
(0.3-30)xl0 5 cm" 3 , CFE of 0.5-60%, non-thermal ve- 
locity dispersions of 0.3-0.7 km s , and rotational-to- 
gravitational energies of 0.003-0.072. For the 18 cores, 
the median number of milllimeter sources within a field 
of view of 0.1 pc is 3.25. 30% of the cores show no 
signs of fragmentation, while ~ 50% split up into > 4 
fragments. In addition, we assessed the number of in- 
frared sources in the field of view, from 2MASS and 
Spitzer data available for the regions. 

- We found no relation between the number of millime- 
ter sources and any of the properties measured for the 
massive dense cores, and especially the correlation seen 
in the optical/infrared between Lboi and the richness of 
a cluster is not seen in our sample of protoclusters. 
Instead, there are hints for a correlation between the 
CFE and the density of the massive dense core, which 
suggests that gravity dominates as density increases. 
Another weak trend is found between the mass of the 
strongest millimeter source and the bolometric lumi- 
nosity, suggesting a 'simultaneous accretion picture', 
from core to fragments, and from fragments to pro- 
tostars. Furthermore, we found a correlation between 
the number of millimeter sources and the number of in- 
frared+millimeter sources, indicating that probably the 
different episodes of star formation keep the same frag- 
mentation level along time. 

- Given the lack of correlation between the fragmen- 
tation level and the core properties in our sample, 
we compared our observations with the radiation 
magneto-hydrodynamic simulations of Commercon et 
al. (2011a) of massive dense cores, and found that 
the low level of fragmentation seen in 122198 or 
120126+4104 can be reproduced with their simulations 
for the highly magnetized core case (with fi = 2), while 
the high level of fragmentation seen in A5142 can be 
reproduced with a dense core with weak magnetization 
(jx = 130). 

Overall, our observations suggest that the optical/infrared 
correlation previously found between the mass of the most 
massive star of the cluster and the cluster richness is not set 
at the early stages of cluster formation but probably arises 
as a result of cluster evolution. In addition, from this study 
we conclude that the fragmentation of massive dense cores, 
rather than being determined by the mass or density of the 
core, could be strongly influenced by its initial magnetic field. 
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APPENDIX 
COMPLEMENTARY MAPS 

As outlined in Sect. 14.11 we compiled a number of millimeter interferometric images published in the literature whose angular 
resolution and sensitivity are good enough to compare to the 4 regions studied in this work with the PdBI in A configuration. 
From these images we estimated the number of millimeter sources within 0.1 pc, N mm . In this Appendix we present the sample 
of 14 millimeter images (see Bontemps et al. 2010 for the millimeter images of the four regions in Cygnus X). 
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Fig. 11. — 1.3 mm continuum high angular resolution maps for cores with 1-2 fragments. Contours for all regions are -4, 4, 8, 16, and 32 times the rms noise 
in each region, listed in TablefJ along with references. Synthesized beams are plotted in the bottom-right corner of each panel, and the black circle corresponds 
to 0.1 pc diameter. Crosses correspond to Spitzer/IRAC sources at 3.6 /im, and plus signs indicate the 2MASS sources. In all panels the red contour corresponds 
to the identification level of 6cr. 
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Fig. 12. — 1.3 mm continuum high angular resolution maps for cores with > 3 fragments (for NGC2071 it is the 3 mm image). Contours for all regions are 
-4, 4, 8, 16, 32, and 64 times the rms noise in each region. The rms noise used for each region is given in Table[3] and references are given in the main text. 
Synthesized beams are plotted in th e bo ttom-right corner of each panel. For 122134, we overplot in green contours the 3.6 cm emission (contours as in Fig.[T}. 
The symbols are the same as in Fig.[TT](red contour corresponds to the identification level of 6cr). We labeled the source used for the SED in cases of possible 
confusion (if no source is labeled, we used the infrared source closest to the strongest millimeter source). 



